地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时，对区内地震在方程中保留震源项，引入正交投影算子进行参数分离，对区外远震采用传统的平滑处理方式，用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理，使得偏导数矩阵中的非零元素成倍增加，对于大型反演问题，这些非零元素常常达到几十GB到几百GB的数量级，巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题，本文研究了偏导数矩阵中非零元素的分布规律，设计出合理的存储结构，采用分布式存储进行矩阵计算，提出了LSQR算法的并行化方案，并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

基于分布式计算模式的GIS的设计与实现

介绍了分布式GIS和分布式计算模式的概念,探讨了基于分布式计算模式的地理信息系统的体系结构、运作机制、关键技术和实现方法,最后对未来的分布式地理信息系统发展进行了展望。     

水文网络模型在分布式流域水文模拟中的应用

水文网络模型是地理信息系统对流域实际河网水系的一种规范化描述,它采用拓扑关系来确定流域水流的空间聚合与分散。准确地描述水文网络对成功模拟流域水流的空间分布起着至关重要的作用。水文网络数据模型的关键在于如何表达水文网络的构成要素,以及如何建立要素间的拓扑关系。本文重点阐述了水文网络模型在分布式流域水文模拟中的应用。实例研究证明,采用水文网络模型来表达流域河网水系结构有助于准确地模拟流域上水流的空间与时间分布。     

基于分布并行的空间数据服务器体系结构分析

从空间数据服务应用需求出发,分析了几种典型高性能存储系统的特点;分析研究了空间数据服务器系统的数据服务过程,提出并设计了由目录服务系统、磁盘服务系统、数据分布系统、客户应用系统以及高速宽带网络组成的分布并行的空间数据服务器系统体系结构。最后,通过空间数据服务实验,验证了该体系结构的可行性。     

A cellular automata model of surface water flow

Previous cellular automata models of surface water flow have been constructed to reflect steady, gradually‐varied flow conditions. While these models are extremely important in showing the near‐equilibrium forms that result from the interactions of water and boundary material, highly dynamic forms and processes require models that represent unsteady flow conditions. In order to simulate unsteady flow in a cellular model of surface water flow, the conservation of mass and the Manning's equations are coupled with an algorithm to delay the movement of water from one pixel to the next until the correct timing is reached. This approach yields highly realistic flood wave hydrographs when compared with flood observations in the Walnut Gulch Experiment Watershed. Coupling this unsteady flow model with simple laws of sediment erosion, transport, and deposition should allow event‐based simulations of watershed and river channel geomorphologic change. Copyright © 2007 John Wiley & Sons, Ltd.     

Airborne gravimetry used in precise geoid computations by ring integration

Two detailed geoids have been computed in the region of North Jutland. The first computation used marine data in the offshore areas. For the second computation the marine data set was replaced by the sparser airborne gravity data resulting from the AGMASCO campaign of September 1996. The results of comparisons of the geoid heights at on-shore geometric control showed that the geoid heights computed from the airborne gravity data matched in precision those computed using the marine data, supporting the view that airborne techniques have enormous potential for mapping those unsurveyed areas between the land-based data and the off-shore marine or altimetrically derived data. Received: 7 July 1997 / Accepted: 22 April 1998     

Taylor定理余项中θ的渐近性和拓广

给出带有m-1项断裂项n+m阶Taylor定理余项中θ的渐近估计,其估计式为 .利用该估计式将Taylor定理推广为 ,其中余项 .最后举例说明该推广在近似计算中的优越性。     

Understanding of hydraulic properties from configurations of stochastically distributed fracture networks

A systematic investigation of the effect of configurations of stochastically distributed fracture networks on hydraulic behaviour for fractured rock masses could provide either quantitative or qualitative correlation between the structural configuration of the fracture network and its corresponding hydraulic behaviour, and enhance our understanding of appropriate application of groundwater flow and contaminant transport modelling in fractured rock masses. In this study, the effect of block sizes, intersection angles of fracture sets, standard deviations of fracture orientation, and fracture densities on directional block hydraulic conductivity and representative elementary volume is systematically investigated in two dimensions by implementing a numerical discrete fracture fluid flow model and incorporating stochastically distributed fracture configurations. It is shown from this investigation that the configuration of a stochastically distributed fracture network has a significant quantitative or qualitative effect on the hydraulic behaviour of fractured rock masses. Compared with the deterministic fracture configurations that have been extensively dealt with in a previous study, this investigation is expected to be more practical and adequate, since fracture geometry parameters are inherently stochastically distributed in the field. Moreover, the methodology and approach presented in this study may be generally applied to any fracture system in investigating the hydraulic behaviours from configurations of the fracture system while establishing a ‘bridge’ from the discrete fracture network flow modelling to equivalent continuum modelling in fractured rock masses. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of linear and nonlinear techniques in river flow forecasting in the Kilombero River basin,Tanzania / Application de techniques linéaires et non-linéaires à la prévision des débits dans le bassin de la Rivière Kilombero,en Tanzanie

Abstract

Modelling of the rainfall–runoff transformation process and routing of river flows in the Kilombero River basin and its five sub-catchments within the Rufiji River basin in Tanzania was undertaken using three system (black-box) models—a simple linear model, a linear perturbation model and a linear varying gain factor model—in their linear transfer function forms. A lumped conceptual model—the soil moisture accounting and routing model—was also applied to the sub-catchments and the basin. The HEC-HMS model, which is a distributed model, was applied only to the entire Kilombero River basin. River discharge, rainfall and potential evaporation data were used as inputs to the appropriate models and it was observed that sometimes the system models performed better than complex hydrological models, especially in large catchments, illustrating the usefulness of using simple black-box models in datascarce situations.